HEPHY-10702 



CP violating asymmetries in single top quark production at the 

Tevatron pp collider 



Ekaterina ChristovaQ 

(N ' 

O ■ Institute of Nuclear Research and Nuclear Energy, 

^ ! Tzarigradsko Chaussee 72, Sofia 1784, Bulgaria 

Oh; 

' Stephan Fichtingerf], Sabine Kraml^, Walter Majerotto^ 

a^ ■ 

Institut fiir Hochenergiephysik der Osterreichischen Akademie der Wissenschaflen, 
>■ \ A-1050 Vienna, Austria 

o 
oo 
O 

Abstract 

o 

:a 

• i-quark production in pp — > W* tb ^ bbW and of the leptons from the subsequent 
(D ! 

r~| , decay W ^ Iv are obtained in the laboratory system. CP violation in the t-production 



X 
J3 



Analytic expressions for the angular distributions of the 6-quarks associated with single 
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1 Introduction 



The experiments at the Tevatron pp colhder offer the possibihty of thoroughly studying the 
properties of the top quark. In particular, in top physics the Standard Model (SM) gives 
negligible CP violating effects [|I| due to the GIM mechanism. Looking for CP violation in top 
quark production and/or decays is therefore one of the best ways to probe New Physics. In 
extensions of the SM as, for instance, in two Higgs doublet models or in supersymmetric (SUSY) 
models, CP violating phases appear rather naturally. They can then cause CP violating effects 
in processes with top quarks at one-loop level. 

In this paper we study CP violating asymmetries in tb and ib production in pp collisions, 
induced by CP violating form factors in the tbW vertex. More precisely, we consider the two 
CP conjugate processes 

pp tb and pp — > — > bt (1) 

followed by the decays t bW^ —>■ bl^ui, i bW~ bl'vi. The possibility of testing CP 
violation in pp tb ^ W^bb was already considered in ^ and in more details in a review 
on this subject in [Q. In this article, we present new additional CP violating asymmetries. In 
particular, we define two types of asymmetries: (i) with b quarks and (ii) with leptons in the 
final states. For the asymmetries with b quarks we assume that the 6's from production can 
be distinguished from those from t decays. This is justified by the requirement that the whole 
event be reconstructed in order to distinguish the events of single t production in from 
the background events . We thus define separate asymmetries for b quarks from production, 
from decay, and from both production and decay. The asymmetries with leptons are in general 
easier to observe. However, the cross section is smaller by the branching ratio W —>■ lu, which 
is roughly 1/3 if all leptons I = e, are counted. 

For both b quarks and leptons, we consider asymmetries in the total number of b and b (/+ and 
/~) and asymmetries in their angular distributions (forward-backward asymmetries). When 
considering the angular distributions of the decay products we assume CP violation only in the 
production vertex of the t quark. As shown in CP violation in t-decay is suppressed by the 



1 



polarization of the t-quark. Moreover, CP violation in the decay would mean that in addition 
to CP violating phases there are also new decay modes of the top quark, which -according to 
the present experimental limits- does not seem to be the case. However, we will also discuss 
the contribution of CP violation in the t decay to the total cross section asymmetries defined 
in Sections 4.2 and 5. 

The paper is organized as follows: In Section ^ we briefly give our notation. In Section |^ we 
define the CP violating asymmetries for b quarks from production, namely total cross section 
asymmetries and forward-backward asymmetries. In Section ^ we derive the formulae for 
the asymmetries for b quarks and leptons from the t decay. We again consider total cross 
section asymmetries and forward-backward asymmetries. The forward-backward asymmetry 
for b quarks from both production and decay is given in Section ^ In Section ^ we give a 
numerical analysis and discussion of the resulting asymmetries in the framework of the Minimal 
Supersymmetric Standard Model (MSSM) with complex parameters. The necessary formulae 
for the form factors, mass matrices, and couplings in the MSSM are given in the Appendices. 

2 Notation 

The quark subprocesses which we consider for single t and i production are 

u + d ^ W ^ t + b, (2) 
d + u ^ W ^ b + i. (3) 

The matrix elements including CP violation read (m^ = 0): 

2 

M'' = ^u{-p;i)^^PLu{pu)-^^u{pt)T''u{-pi), (4) 
2 s- 

M"' = ZE_u{-pu)laPLu{pd) . ' , n(pfc)r"^(-p,-) (5) 
2 s - mfy 

with 

P" = ^-p^(l + ^fCP) + £lp^^gCP^ (6) 

nit 



= — Pj, -P" = Pb ~ P? -Pl = (1 — 75)/2. /f^''' and 5'^''' are the CP violating form 
factors of the tbW vertex. They are complex functions. The asymmetries considered in this 
paper measure the absorptive parts of these form factors — f^^ and ^mg'^^. The real 
part of (7^^, ^eg'^^, can be measured through triple product correlations. As there is no CP 
violation at tree level, 3fJe f^^ has no physical meaning and cannot enter measurable quantities. 



3 b quarks from production 

In the centre-of-mass system (CMS) of u and d 01 u and d, with the z axis pointing along the 

momentum pu or pdi the angular distributions of the h and h quarks in the subprocesses @) 
and (0) are: 

dcrf = C[a^o + a'l cos 61 + a\ cos^ 61 } d cos 61 , (8) 

daf = C[al-a\ cos 61 + a\ cos^ 61 } d cos . (9) 

Here 

a? = ar±ar, (10) 

a'o'' = -K + 5)/2, (11) 



,5M 



and 



(12) 



af^ = {m',-s)/2, (13) 
a^^ = 2af^53m/f^(S) + (S-m,2)!3m^?g^(S), (14) 

^CP _ o ^SM /CP 



2ar53m/r(5), (15) 
2af ^ (53m /r(5) + 53m ^^^^(S)) , (16) 



16[s(s-m2,)]2 ■ ^ ^ 

s stands for 

S- = iPu+PdY for a*^ 

(18) 

s+ = (p„+pj)^ for a*'' 



The polar angle cos 9 in the laboratory frame — the CMS of p and p — is related to the angle 
cos 9* in the CMS of the initial quarks by: 

cos9 -V 
1 — V cos 9 

where v is the velocity of the laboratory system in 

Here cos^ stands for cos^;, or cos9i, and Xi and X2 are the '^^^^^.^^^^ ^^^^ ^^^^f^^^^^.^^^. ^^^^^^^^^^^^ 

' T . -1 . • r ^j^^ ^ ^ quarks in the laboratory system thus is 



of the quarks. The angular distribution of 
^ = d£2^. We hence obtain: 

a cos a a cosy* rfcosS 



(19) 

ithe CMS of the quarks, v = (xi— X2)/(xi+a;2). 
^2 are the fractions of the longitudinal momenta 



.tb 



da' 
d cos 9ii 



{l-vcos9bY 



daf 



.tb 



C{1 



dcos9i (1 — f cos6'fe)' 



00(1 — V cos 9i,y 

+ a\{cos9b — v){l — vcos9h) + 03 (cos 6*;, — v) 
ao(l — V cos 6^^)^ 

— a5(cos6'5 — f)(l — wcos6'5) + a2{cos9i — v) 



(20) 



(21) 



Taking into account the two possibilities xi = x„, X2 = Xd and xi = Xd, X2 = Xu, and using the 
parton distribution functions /„ and fd of the proton, we get in the laboratory system: 

daf 1 fdaf{s,v) 



d cos 9b 

daf 
d cos 9-h 



2 J d cos 9b 
1 fdaf{s,v) 



2 J d cos 9u 



[fu{Xu)fd{Xd) + fu{Xd)fd{Xu)] dXudXd, 
[fu{Xu)fd{Xd) + fu{Xd)fd{Xu)\ dXudXd 



(22) 
(23) 



with 



s = XuXdS and v = {xu- Xd)/{xu + Xd) 



(24) 



3.1 Total cross section asymmetries for h from production 

For the total cross sections of (0) and (^ we obtain: 
{a,f'' = 2C(a^'"^ + af/3) 



SM/ 



CP _L 



1 ± 2 /^^ ± 



2s 



(25) 



where 



m 



2\2i 



ml + 2s) 



(26) 



2A[s{s-ml,)f 

and aw = g^j^'R. In (^) the upper sign stands for and the lower one for a*. The CP 
violating asymmetry for the total number of h and h quarks from the production processes ([l|) 
is thus given by 



From fE2 



nr{b) 



3f) and (^) we obtain: 

/ C (a'^^ + a^^/3j fu{xu)fd{xd) dxu dxd 
J C (a^^^ + af 73) Uxu)Uxd) dx., dx^ 



(27) 



2/a 



SMi 



^mf^^{s) + ^^^f^$5m5()J^(s)J fuixu)fdixd) dxudx 



(28) 



/ d^^'{s) fuixu)fdixd) dXu dxd 
This asymmetry was already obtained and discussed in Q]. We recall it here for completeness 
and also include it in our numerical analysis. 



3.2 Forward— backward asymmetries for b from production 

We next define two CP violating forward-backward asymmetries for the production process: 



(^i)f + (^i)'s 



and 



I - {^1)3] [fu{Xu)fd{Xd) + fu{Xd)fd{Xu)] dXudXd 

I li^lfp + (o'i)b] [fu{Xu)fd{Xd) + fu{Xd)fd{Xu)\ dXudXd 



(^1)^ + (^i)'f 



/ [{^i)b - {^i)f] [fu{Xu)fd{Xd) + fu{Xd)fd{Xu)] dxudxa 



(29) 



(30) 



/ [{^l?B + (^i)f] [fu{Xu)fd{Xd) + fu{Xd)fd{Xu)\ dXudXd 

Here {(TiYf{b) number of b quarks from production in the forward (backward) direction. 

Note that in single t production, contrary to ti pair production, we have \Af^\ ^ l^f^'^l- This 



5 



is due to the fact that CP violation leads to a difference in the total cross sections of b and b 
production. In terms of form factors we obtain: 



A^ib) 



with 

h{s) 



J C (a^P + (1 - v^) af ^/2 + a^^/s) fu{xu)fd{xd) dx^ dxd 
J C (af^ + (1 - v^) aPV2 + af 73) dx^ dxd 

2 / h{s) [(2m2 + s{7 - ?,v^)) /f ^ + 2(m? - s) S^m g^^] /„(x„)/,(a;,) dx^ dx^ 



f h{s) [2mj + s(7 - 3t;2)] fu{xu)fd{xd) dxu dxd 

JC (a^^ - (1 - v^) af^/2 + a^^/s) fu{xu)fd{xd) dx^ dxd 
J C {ai^' - (1 - v^) af*72 + af 73) fu{xu)fd{xd) dx^ dxd 



(31) 



2 / /z(s) [(2m2 + si3v^ + 1)) 53m /^^ + 2(m2 - S) g^^] /„(x„)/d(x,) dx„ 



/ h{s) [2ml + ■s(3t;2 + 1)] fu{xu)fd{xd) dxu dxa 



s{s - m^) 



(32) 



(33) 



4 Secondary b quarks and leptons 

Let us now turn to the b quarks and leptons which originate from the t decays: 

t-^bW, bl+u . (34) 

Following [0, in the narrow width approximation {Ft rrit) we obtain in the CMS of ud: 



d&n 



' da\ \ dV{t 
— a cos Of 



X. 



tot 



El 



'^''*-) d cose; "^{t ^ X + ...) Br{t ^ xX) 



dcose*t J ' r 



(35) 



and analogously for the decay of t. In (|35D, dcr{/d cos 6^ is the distribution of the t quarks in 



dd\ 
d cos 61 



C[al + a\ cose; + a\ cos^f } , 



(36) 



and E* is the energy of the decaying t. ^(t —^x + ...) is the angular distribution of the 
secondary particle x {x = b,l~^) in this frame with t polarized, normalized to the partial decay 
width r(t — s> x). The branching ratio Br{t — > xX) stands for Br{t bW'^) or Br{t bl^ui). 
If /3j is the velocity of t in the CMS of ud we have 0: 

V r ; 4vr Ef (1 - A* cos Bt:f \ (ptP.) j ' ^ ^ 

where determines the sensitivity of the particle x = b,l to the polarization of the t quark: 



ai = -1. (39) 



^ is the polarization four-vector of t, which is determined by the production process and 
cos6'^2, is the angle between the momenta pt and Px'. 

cos 6*^ = sin e* sin 6** cos 0* + cos 9* cos 6*; . (40) 

The treatment of the t polarization four-vector and the general formula for the differential cross 
section, in the CMS of (ud), in terms of the SM- and CP-violating components of ^, needed to 
derive the angular distribution daf / d cos 9x are given in Appendix 0. 

4.1 Angular distribution of the decay products 

Integrating (^) over dcos9^ and dcpl we obtain the angular cos^* distribution of the decay 
products X = b,l~^ in the CMS of ud: 
dcT^ 

i3[6o + &icos^: + 62Cos2C] (41) 



d cos 9% 
where 
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For the CP conserving part we obtain: 



SM 
^2 



"o 



SM 



dl"" 



1 



2m1 



mt 



2mt 



mt 



2sm^ [m^ + s) In — (m^ + 3sm^ — 5s + s" 



-2sm; In — ^ 

s 



4 -2 

— s 



— 6m^ ^(mj + s) In — — h (m^ — s) {rri 



[ml + s) In ^ - 2{ml - s) 



2s 



-3d: 



s in — (m^ — s) 

s ^ 

SM 



lOsml + s^) 



(45) 

(46) 

(47) 

(48) 

(49) 
(50) 



For tlie CP violating part we have: 



„CP 



c7 



^CP 



d\ 



'CP 



dr 



^SM fCP , "^f + 



9 , "ir 
-2s In — 



+ — 



CP 

ksmg)^ . 



ml 



mr 



-6m^ s (m^ + s) In — h + 9s m^ (m^ — s) — s 

5 



CP 



+ 



2 ml {ml + 2s) In 



5mf + As ml + s^ 



2s 



-2rfr53m/r + ^ 



—2s m, In h m. — s 

s 



o CP 



o CP 
3m , 



(51) 
(52) 

(53) 

(54) 

(55) 
(56) 



The coefficient B is: 

^ —nal^ml Br{t ^ xX) 



(57) 



8[s(s-m2^,)]2 • 

Relations (|50| ) and (|56D ensure that the polarization of the t does not contribute to the total 
number of b quarks from t decay. The angular distribution in the laboratory system is obtained 
from (^iD by the Lorentz boost (|I9D: 



d cos 9r 



1 r do^'^ (s v) 

77 / —r [fuiXu)fdiXd) + fu{Xd)fd{Xu)] dXudXa 

z J a cos Or 



(5J 



with 



dcosOx (1 — f cos6'2.)'^ 



60(1 -vcosdr,)'^ 

+ bi{cos9x — v){l — V cosOx) +b2{cos9x — vY , (59) 



and s and v defined in (|2^) . The angular distribution of the decay products from t decay is 



obtained from (|59D by CP conjugation. 

4.2 Total cross section asymmetries for the decay products 

Analogously to ( p8D we define a total cross section asymmetry TZ^2^{x) for the secondary particles 
x: 

7^f(a;) = ^^Z^, x = bj+. (60) 

This asymmetry was first suggested for b quarks in [Q, but without giving an explicit analytic 
expression. Using our result (|^) and ( |59D we obtain: 

/ B{b'^^ + b^^ /'i)fu{xu)fd{xd) dxu dxd 



' ^ ^ J B{bi^^ + 6f^/3) A(xJ/,(xrf) dx^ dxd 

J B{c^^ + c^^/3)fu{xu)fd{xd) dxu dxd . 

(61) 



/ B{cl'' + cf ■V3)/.(x„)/rf(x,) dx^ dxd ■ 
As there is no dependence on the t polarization, ( |6lD implies that the total cross section 

asymmetries for secondary b quarks and leptons are equal in magnitude: 

nrib) = K\i). (62) 

This is a consequence of CP invariance in the decay W —>■ Iv. Eq. (^) also implies that 

7^f (x) = -7^f (6), X = 6, /+ . (63) 

However, ( |63|) is not valid in general. As already pointed out in 7^2°* (^) receives contribu- 
tions from both the production and decay vertices. In the general case, (|63|) reads: 

7^f (x) = -7^f (6) + (5f ^, x = b,l+. (64) 



where 6^'^ is due to CP violation in the t-decay fl^ 



CP _ T{t bW+) - T{t bW- 



^ Y{t^bW^)^Y{t^bW~) ^ ' 

Testing (p^ would be a model independent test of CP violation in the t-decay vertex. 
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4.3 Forward— backward asymmetries for the secondary products 

For the secondary products x we can again define two forward-backward asymmetries in tlie 
laboratory frame: 

{(72)f - ((^2)1 / [(0-2)^ - {^2)%] fu{Xu)fd{Xd) dxu dXd 



R^ix) 



R 



FB, 



X] 



From (|T) and 



witfi 



his) 



his) 



and 



R^ix) 



witli 



(0-2)f + {(^2)% I [i^2)F + {^2)3] fu{Xu)fd{Xd) dXu dXd ' 

(q'2)| - (q"2)| ^ / [((72)1 - (0-2)!] fu{Xu)fd{Xd) dXu dxg 
(0-2)1 + (0-2)1 / [((5-2)| + ((5-2)1] fu{Xu)fd{Xd) dXu dXd ' 

I) we obtain: 

/ B (c^P + c^P/3 - {v^ - 1)6^/2) Ux^)Uxd) dXu dx^ 



J 8 (co^*^ + cf^VS - {v^ - l)6f 72) fu{xu)fd{xd) dxu dxd 

-2jD{s) {ki{s) ^mfi + h{s) Sm^a} fu{xu)fd{xd) dXu dxd 
J D{s) ki{s) fu{xu)fd{xd) dXu dxd 



{ml - s) 2n4 3t;^) - sm^((l - 3t;^) - 6 a^.(l - v'^)) 



+ a^){l-v^)mi r In 



[m] - s) 2{ml - s)^ + 3 a^s{l - v"^) [m] + s) 



— 6 ax{l — V )m^s In 



J B (djP + c^^/3 + (t;^ - IW/2) Uxu)fd{xd) dXu dxd 
J B icl'^ + 4^IZ + {v^ - l)6f 72) Uxu)fd{xd) dxu dxd 

-2 j D{s){k2{s) fi + /2(g) '^mgR}fu{xu)fd{xd) dXu dxd 
JD{s) k2{s) fu{xu)fd{xd) dXu dxd 



2mt - s2(l + 3v^) + sml{{5 - 3v^) - 6a^(l - v^)) 

2 

mr 



2{ml - sf - 3a^.(l - v'^) s {ml + s] 



+ 6 axil — V s In 



mr 



(66) 
(67) 



(6J 



(69) 



(70) 



(71) 



(72) 



(73) 
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The factor D{s) is 

D{s) = [s{s - m'^r' 



(74) 



Note that i?f^ are polarization asymmetries, i.e. they measure different combinations of the 
CP violating contributions to the t quark polarization. This can also be seen from the explicit 
expressions for ki^2 and li^2- As a consequence, the forward-backward asymmetries for the b 
quarks are different from those for the leptons. Note, moreover, that the contribution from 
t polarization enters only through the term bi which is linear in cos^-^, which ensures that 
polarization does not contribute to the total cross section. 



5 b quark forward— backward asymmetry for the sum of 
the cross sections 

An asymmetry that seems most convenient what concerns statistics is the forward-backward 
asymmetry for the total number of b and b quarks from both production and decay: 



■rb _i_ J) 



(7°p + (7? 



where 



6,6 / \o.o I / \o.o 0.0 / \o.o I / \o,o 

(^F = + K(^2)f : CTfi = + {'^2)b ■ 



\6,6 



\b,b 



b,b 



\b,b 



\b,b 



We have: 
A^'^ib) 



(75) 



(76) 



(77) 



/ [4 C (ao^^ + + (1 - v^) {C af^ + B 6D] /„(x„)/,(x,) dx^ dx, ' 

where we have used the fact that in the SM the total numbers of b quarks from production and 
decay are the same: 



(C a^^ + B 6^^) + (C af ^ + B b^^) /3 = 2C(a^^ + af ^/3) . 

Explicitly, in terms of the CP violating form factors, wc obtain: 

pj^^^^ QJD{s)\k{s)'^mf2^ + l{s)'^mg'^^\ fu{xu)fd{xd) dxudxa 



A'^'ib) 



J D{s) m{s) fu{xu) fd{xd)dxu dxa 
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(78) 



(79) 



with 



k{s) 

l{s) 
m(s) 



2s 



(s — mj)(s + (2 — Sabjm^) + 3(1 — a^jm^s In 



m; 



-2 



o 2-1 

2mj s m 



+ 3m,V(l-a,)(l-^;2)ln^ 



(80) 
(81) 

(82) 



Assuming CP violation only in the t-production vertex, the total cross sections cr* and a'^ 
for the b and b quarks from both production and decay are equal, i.e. the total cross section 
asymmetry 7?.*°* (6) is zero: 



a' -a" 



= • 



However, when CP violation in t-decay is also considered, we have: 

sjCP 

TV°\b) = \- ■ 



(83) 



(84) 



In the general case, we thus have three different total cross section asymmetries and they mea- 
sure three different quantities. The total cross section asymmetry of the 6-quarks from produc- 
tion, 7^^°* (6), measures CP violation in the production vertex, the total cross section asymmetry 
of the i-decay products, TZ^2*{x), x — b,l measures CP violation both in the production and 
decay vertices, the total cross section asymmetry of the 6-quarks from both production and 
decay, 7?.*°* (6), measures CP violation in the i-decay vertex. 

Neglecting CP violation in f-decay, we are left with only one forward-backward asymmetry 
A^^ when all b and b quarks from single t-quark production are counted. 



6 Numerical results and discussion within the MSSM 

In this Section we present numerical results for the discussed asymmetries in the framework 
of the Minimal Supersymmetric Standard Model (MSSM) with CP violating phases. In the 
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MSSM one has the gaugino mass parameters Mi = M', M2 = M, M3 = rrig corresponding to 
the groups U(l), SU(2), SU(3), where nig is the gluino mass. We assume the GUT relations: 



(a,/a2)M~3M, M' = (5/3) tan^ ^v^M ~ | M (85) 



with a common phase of Mj which can be rotated away by an R rotation. The only complex 
parameters relevant for our discussion thus are: the Higgsino mass parameter /i = e*'^'' and 
the SUSY breaking trilinear couplings of the stops and sbottoms, At = \At\e^'^* and = 
\Ai) \ e^'^''. As we work in the limit of = the phase (pb does not play any role. The phase of /i 
is strongly constrained by the upper bounds of the electric dipole moments of the electron and 
neutron. 0^ must be very small except the SUSY masses are large (> 1 TeV) or there are 



strong cancellations between the different contributions |]I2[. We therefore take 0^ = 0. After 
this the only relevant phase we are left with for the considered CP violating asymmetries is the 
phase (pf 

In the numerical analysis we assume CP violation only in the production vertex and not 
in the decay of the t quark: CP violation in the decay would mean that in addition to CP 
violating phases there are also new decay modes of the top quark. According to the present 
experimental limits there is only a small window left for the decay t Xi ~^ ii- A detailed 
analysis of CP violation in t-decay was performed in 

At one-loop level, the reactions ub —>■ tb, du — > bi receive radiative corrections from triangle 
and box graphs with charginos xt ■> neutralinos x% squarks qi {i = 1,2; j = 1...4), and gluinos 
g in the loops. The analytic expressions for the form factors due to these diagrams have been 
worked out in [Q. Following [Q, in the limit 0^ — 0, rriu^d — ^ 0, only the graphs with (x'^x'^t) 
and {t b g) loops (see Fig. |l]) contribute to CP violation and the contribution from the box 
diagrams is negligible. We thus base our numerical analysis on the two contributions - from 
the {x^x^t) and {tbg) loops. Our formulae for the form factors agree with those in P] and are 
given in Appendix p. 

We now want to analyze the influence of a possibly large phase of At on the CP violating 
asymmetries defined in the previous Sections. For this purpose we choose three scenarios 
of gaugino-higgsino mixing: (i) a gaugino scenario with M = 116 GeV, /i = 400 GeV, (ii) 
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a Higgsino scenario with M = 400 GeV, fi = 116 GeV, and (iii) a 'mixed' scenario with 
M = fi = 168 GeV. Further, we take tan/3 = 4, = 300 GeV, = 270 GeV, and 
mi) = 330 GeV. The SM parameters are: rrit = 174 GeV, m^y = 80.03 GeV, sin^ 9\y = 0.23 
a{mz) = 1/129, and ^^(m^) = 0.12 (we neglect the bottom mass, = 0). For the structure 
functions of the proton we use the CTEQ5 parton distribution function cteqSml [|T3] . We leave 
out the theoretical uncertainties associated with the QCD corrections and parton distribution 



functions discussed in |T^ since these uncertainties, being CP-even, cannot mimic the CP- 
violating asymmetry discussed here. 

Figures H-^jshow the resulting CP violating asymmetries for the three scenarii as a function 
of (f)f \At\ is chosen such that mi ^ 96 GeV. Moreover, we have m-+ ~ 104 GeV and m^o ~ 50- 

1 Al '^l 

100 GeV. The asymmetries are typically of the order of 10~^. Note, however, that the masses 
in our scenarii are just at the borderline of the experimentally allowed values [^. If the mass 



spectrum becomes heavier the CP violating asymmetries quickly decrease. The asymmetries 
also decrease with increasing tan (3. 

The asymmetries A2^{b) and -Rf^(fe) seem to be the most promising ones — they reach up 
to ~ 0.1%. Our results are an order of magnitude smaller then the estimates obtained in 0]. 
A general result for the two types of forward-backward asymmetries is that 

A^^ > nf > A(^, R^^ > i?f«, (86) 

and A2^ {^2^) some cases can be ~ 2 — 3 times bigger than Af^ (Rf^). Another general 
feature is that the forward-backward asymmetries for the decay products i?f(f) are almost equal 
for &-quarks and leptons: 

i?f(f) (b) ^ i?f(§ (/) . (87) 

The forward-backward asymmetry A^^ for the sum of b quarks from both production and decay 
turns out to be extremely small < 0.02%. 
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7 Conclusions 



We have studied the process of single i-quark production in pp coUisions. Assuming CP violation 
in the ^-production vertex we have defined different angular and total-cross-section asymmetries 
as measurable quantities. General analytic expressions for these asymmetries in terms of the 
corresponding CP violating form factors are obtained. Relations sensitive to CP violation in 
the t-decay vertex are defined. We want to emphasize that the formulae are valid independently 
of the origin of CP violation. 

We have performed a numerical analysis within the MSSM with complex phases. In accor- 
dance with the constraints from the measurements of the electric dipole moments of the electron 
and the neutron we have only kept the influence of the phase of Af. The effects turn out to be 
rather small. The discussed asymmetries are of the order 10~^ - 10""^ and they decrease as the 
mass spectrum and tan (3 increase. 

The cross section for pp — > W'^ — > i6 is ~ 340 fb. This means ~ 10^ events for RUN II of 
the Tevatron with an integrated luminosity of 30/6^^ (15 per experiment). To measure 
an asymmetry of the order of 10"*^, one would need about 10^ events. This is far beyond the 
observability of such an asymmetry at the Tevatron. Nevertheless, it would be worthwhile to 
look for CP violation in top physics as it would imply not only physics beyond the Standard 
Model, but also beyond the Minimal Supersymmetric Standard Model. 
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A Differential cross section for the t decay 



A.l The t polarization four— vector 

Let us write the amplitude of the process (^) in the form: 

M = u{pt) Ml u{-pi) 
where 



Ml = ^ w(-pj) 7a u{pu) r° . (89) 

2 2 s — 

The polarization four-vector of the t quarks is determined by 

, ( Pip'A Tr{MiA(-p-,)MiA(p,)7V} , ^ 

V ^? / Tr{MiA(-p5)MiA(j9t)} 

In the most general form the covariant decomposition of reads 

r = Qr^i + Q2H + e{a,p^,p3,pt)V\ (91) 

Here 

Q1 = p:-^^p?, Q2=Ps-^^pt (92) 

are two four-vectors in the production plane orthogonal to pt, and e(a, Pu, Pd, Pt) is a four-vector 
orthogonal to pt and to the production plane. Using the notation 

U{ct) =al + a\ct + a\ , q = cos 9t , (93) 

and the usual Mandelstam variables 

s = {Pu+Pd)^, i= {Pu-Pbf = {Pd-Ptf , s + i+u = mf (94) 
we obtain in the CMS of ud: 

- u{ml -u) + ml{Au + 3s) - } , (95) 
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{-4mtt(l-253m/f^) 



CP 

1 



(m^ - i) \-2{ml - u) - t] - u{m^ - u) - m^{Ai - 3s) 



— s 



[mr 



nit 



CP 



[s — m. 



(96) 
(97) 



A. 2 General formula for the differential cross section 

From (^5|) , (W^) and (|9T|) we obtain the general formula for the differential cross section for the 
production of the t quarks and their subsequent decay in the CMS of ud: 



C_ m^^Brjt xX r , srcp 

47rEKl-Acose,.)2 V 



A-n 
+ a 



rritVs 



2EA1- PtCosOt,] 



(Pf + V^"") ( cos 9, - Pt cos 9t^il - Pt cos Ot. 



1 - 



mi 



[1 - (3t cos 9 tx) 



mt 



where we have used the notation 



d cos 9t d cos 9x dcl)^ 



v{±vl, V± = vi^' + ^ 



Ar(Q) 



V 



met) 



and 



Af{ct)=Af'''{ct)+Af^^{ct) 



Et 



s + mi 



\Pt\ 



s — m 



{nufLtfix) denotes the triple product 
{nuUtfix) = Hu X fit ■ fix 



(3t 



s — 

s + ml 



(98) 

(99) 

(100) 
(101) 



(102) 



with fiu the unit vector pointing in the direction of etc. 
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B The CP violating form factors in the MSSM 



B.l Chargino— neutralino— stop loop 



The CP violating form factors from the t„ xtxi J = 1? 2; /c = 1...4) loop are: 



+ O2 mtrhl C2 - C3 mtm^iCo + C2) - O4 m+m° Co} 



(103) 



and 



nit 



{Ornit Cu - 02ml C2 + OmfiCo + Ci + C2)} 



(104) 



Here we use the notation = and = iTi^o,. The Cx are the standard three-point 
functions for which we follow the convention of [l7|. In this case, 



Cx = Cx{ml ml s, m|+, m? m|o) , X e {0,1,2,00,11,12,22}. 



See Appendix p.3| for the explicit definition of the Cx- The couplings Oi = op-'^ are 
Ol^^' = -V2 -smiling J\ 



nj kj nkJ ' 

njk 



'3 

= -y2 9m(4.0ga*;). 



(105) 

(106) 
(107) 
(108) 
(109) 



with ll^j, Oj^j^, al^j., and 6^;. given in Appendix |C|. Notice that in (|103D and (|104|) one has to 
sum over n,j, k. 



B.2 Gluino— stop— sbottom loop 



The g inbm loop gives: 



/F^(^) = o, 



2a 



9r {g) = --—mtm-gO^Ci 

6 71 



;iio) 
'111) 
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where 

Ci = Ci{ml m\, s, m?^, m?, m?J (112) 

and 

o, ^ or = irLr ^ra (r^ry . (113) 

With the exphcit expression for the stop mixing matrix we have, see Appendix |C|, 

irLl' = {cos2%,sin2^a, m = l,2, (114) 
53m(rtrL) = ^sin^,-sin2e,-{l,-l}, n=l,2. (115) 

Again, in ( |111| ) a summation over m, n is assumed. 



B.3 Three point functions 



Here we give the definition of the Passarino-Veltman three point functions |T^ used above in 



the convention of |17|. For the general denominators we use the notation 



= _ and = {q + pjY - . (116) 

Then the loop integrals in D = 4 — e dimensions are as follows: 

1 



Co = 








ZTT^ J 




c, = 






in^ J 






m"^ J 





(117) 

PimCi+P2mC2, (118) 



2)0plp2 

= g^luCQQ + Pl^^PluCii + {pif,P2u + P2^lPlu)Cl2 + P2f,P2uC22 ■ (1 19) 

where the C's have {pi, {pi — p2y , pi, ml, ml, 1712) as their arguments. 



C Masses, Mixing Matrices, and Couplings 

The neutralino mass matrix in the basis of 

vl/0= (-^A^-^A^V^°,^,V^°,J (120) 
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is: 



M' —mz sin 6^/ cos /9 mz sin ^^vi^ sin /? 

M mz cos 6*^ cos (3 —mz cos 6*1^ sin /3 

— mz sin ^vt^ cos (3 mz cos cos /3 — 

mz sin sin /3 —mz cos sin /3 -/i 

Here tan/5 = V2/V1. This matrix is diagonalized by the unitary neutrahno mixing matrix N: 



121) 



where Al^ is a diagonal matrix with non-negative elements 
of the physical neutralino states. 
The chargino mass matrix is: 

M V2mwsm(3 ^ 

y/2mw cos P n 
It is diagonalized by the two unitary matrices U and V: 



rhi, 777.2, "^3) "^4 — masses 



;i22) 



where is a diagonal matrix with non-negative entries — 777.1", 777.2" — the masses of the 
physical chargino states. 

The mass matrix of the stops in the basis (^l, in) is 



I + "^|cos2/?(i - |sin^^iy) -F7?7f 



{At - /X cot /5) 777* 



(^4* — ij, cot P) mt 
777^ + ^m% COS 2(3 sin^ 6'^' + 

Mf is diagonahzed by the rotation matrix F* such that F*^ Al~ F* = diag{m'~^, m~^) and 
F*(|). We have: 



(123) 



F* = 



F* F* 



( e-^ico^di -e-^'^*-sin^j- ^ 



ea"^* sin^f 



6 2"^* COS^^ 



(124) 
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The interaction Lagrangian which we need is: 

Af„x; = 9b{ki^PL+lljPR){xtrL + h.c., (125) 

A?„^« = ^t(6LPL + <,PR)x^~n + h.c., (126) 

^xtxlw = 9 xl r {O^.Pl + Og.P^) x; + h.c. . (127) 



The chargino-stop-bottom couphngs are: 

4- = -^.irL + ^*^.2rL , (129) 

with 

v2m^ysinp v2mvKCOsp 
The neutrahno-stop-top couphngs are: 

^nk ~ /Ife + ^Lfc Tfin ' (131) 

Kik — Tin + Ink Tiin > (132) 

with 

fik = -7i(iVfc2 + |tan^^iVfci), (133) 

/J, - ^tan^H^TV*,, (134) 

/iL = -^tAT*, (135) 

in the basis = (— ^A', — «A^, -0^^, -0^2)- The chargino-neutrahno-H^ couphngs are given by: 

Okj = -^^NkAV;^ + Nk2V;,, (136) 
= ^^:3t^.-2 + A^:2t^.-i- (137) 
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Figure 1: One-loop Feynman diagrams for ud tb. 
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Figure 2: CP violating asymmetries (in units of 10 ^) for the gaugino scenario: M — 116 GeV, 
^l = 400 GeV, tan/3 = 4, = 300 GeV, rrifj = 270 GeV, = 330 GeV, ~ 104 GeV, 
and rrif, ~ 96 GeV. 
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Figure 3: CP violating asymmetries (in units of 10 ^) for the higgsino scenario: M — 400 GeV, 
^l = 116 GeV, tan/3 = 4, = 300 GeV, = 270 GeV, = 330 GeV, m-+ ~ 104 GeV, 



and mi^ ~ 96 GeV. 
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Figure 4: CP violating asymmetries (in units of 10 ^) for the 'mixed' scenario: M — jj, — 
168 GeV, tan/3 = 4, = 300 GeV, = 270 GeV, = 330 GeV, m-+ ~ 104 GeV, and 
mf, ~ 96 GeV. 
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